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Complete spectral gap in coupled dielectric waveguides embedded into metal 
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We study a plasmonic coupler involving backward (TMqi) and forward (HEn) modes of dielectric waveguides 
embedded into infinite metal. The simultaneously achievable contradirectional energy flows and codirectional 
wavevectors in different channels lead to a spectral gap, despite the absence of periodic structures along the 
waveguide. We demonstrate that a complete spectral gap can be achieved in a symmetric structure composed 
of four coupled waveguides. 



Negative index matamaterials (NIM) are artificial ma- 
terials which have simultaneously negative permittivity 
and negative permeability.^^ In NIM waveguide, modes 
are backward when more energy flows in NIM than in 
other channels. Coupling of a forward propagating mode 
in a conventional dielectric waveguide with a backward 
mode in a NIM waveguide has been investigated theo- 
retically in both linear and nonlinear regimes. ^^^ When 
a forward mode is coupled to a backward one, the back- 
ward mode transports energy backwards, leading to the 
formation of spectral gaps without periodic structures 
along the waveguide. This feedback mechanism may play 
an important role in nanophotonics, as it could signifi- 
cantly simplify complex geometries that are required for 
subwavelength optical manipulation and concentration. 
However, due to the fabrication complexity and high 
losses of NIM, coupling involving NIM is currently not 
experimentally feasible and therefore the mechanism has 
not attracted significant attention. 

There has been a surging interest in the field of plas- 
monics, as it offers one of the most promising approaches 
for subwavelength optical concentration and manipula- 
tion (for a comprehensive review, see e.g. Refs^— ). 
In some plasmonic structures, backward modes exist in 
regimes when more energy flows in the metal than in the 
dielectric. ^^^^. These structures are much simpler and 
more fabricable than those involving NIM. In this let- 
ter, we propose a design of plasmonic coupler involving 
the coupling between the backward TMqi and the for- 
ward HEii modes in dielectric waveguides embedded into 
metal [see Fig. [T]^a)]. We find a polarization dependent 
spectral gap in a structure of two coupled waveguides 
and a complete polarization independent gap in a Cs^ 
structure with four coupled waveguides. 

It was recently reported^^ that taking experimental 
data of bulk metali^ in numerical calculations of plas- 
monic modes may lead to losses which are much higher 
than real losses observed in experiments. In our study, 
we use the Drude model to simulate the optical proper- 
ties of a metal: £m{^) = 1 — oJp/oj{uj + iujr)^ where ujp is 



the plasma frequency and oOr is the collision frequency. 
At the same time, we define two normalized quantities: 
loss 7 = ujr/ojp^ and size parameter a = Rujp/c^ where R 
is the radius of the dielectric core, and c is speed of light. 
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FIG. 1. (Color online) (a) Two dielectric waveduides with 
£i = 9 and S2 = 4 separated by D embedded into infinite 
metal. Green arrows indicate the energy flow at different 
channels for the wavevector along z; (b) dispersion of two 
coupled waveguides. Yellow region indicates the incomplete 
polarization dependent spectral gap obtained using temporal 
coupled-mode theory with Vg^O.lSc, Vp3=0.039c, and ^=0. 



Fig. [TJa) shows the two- waveguide structure we study: 
two dielectric rods of the same radius a = 1.21 (corre- 
sponding R is about 25nm for silver) but with £1=9 
and £2 = 4 embedded into infinite metal. First, by 
analysing the dispersion of a single waveguide, we find 
that the backward TMqi mode for £ = 9 intersects with 
the forward HEn mode for e = 4 at uj/ujp = 0.3856 [see 
FiglSfa)]. This point corresponds to A ^ 400nm for sil- 
ver. For the TMqi mode, more energy flows in the metal 
than in the dielectric, which is similar to the backward 
SPP on metallic wires ^i^ It should be emphasized that 
the directionality of TMqi and HEn modes are radius de- 
pendent: the TMqi mode can become forward when the 
radius increases, and the HEn mode can become back- 
ward when the radius decreases. ^^ However, the HEn 
mode has linear polarization inside the dielectric [see 
Fig. EKc)], which could be excited directly with a nor- 
mal incident wave^^, whereas the TMqi mode has radial 
polarization [see Fig.[2fd)] with much higher losses in the 
coupling region [see Fig.[2fb)]. Prior to numerical study, 
we use temporal coupled-mode theor y -^^i^^ (TCMT) to 
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FIG. 2. (Color online) (a) Dispersion curve and (b) losses of 
TMoi mode for £d = 9 and HEn mode for Sd = 4. (c) and (d) 
Poynting vector component Sz (colourmap) and transverse 
electric field Et (arrows) for HEn mode and TMoi respec- 
tively at a;/cjp=0.3856 with 7=0.002. 



get a qualitative understanding of dispersion relation in 
the lossless case. The eigenmodes of a coupled system 
are expressed as a superposition of individual waveguide 
modes: E = ^^ A^(z, t)E^(x, ?/)e^'^'^""0+'^'^'^^^ where 
^mo = ^m at /c = /co and Hirnm IS sclf coupliug Coefficient. 
For the two coupled waveguides, three modes can couple 
to one another: two forward HEn modes of preferred x 
and y polarizations, which could be approximately recon- 
structed by two orthogonal eigenmodes of circular polar- 
izations: Ai^2{z,t), and backward TMoi mode: As{z^t). 
The coupled- mode equations in time domain are: 

. ^^\^^^-^^ +^.. ^^^-f^^-^^ +.A3(.,t)e-^^^=0 



dt 



dz 



.dAs{z,t) , . dAs{z,t) 



dt 



IVgS- 



dz 



Y,Am{z,t)e-'^''=0. 



where S = |(a>:33 + cuso - ^ii - ^lo) = |(/^33 + ^30 - 
/^22 — ^20 ) is the antisymmetry parameter of two waveg- 
uides; ^1,2,3 are normalized envelopes; Vgi {vg = '^^1,^2 > 
0, Vgs < 0) are the group velocities at ujq = uj{ko); /t:i2 
= A>:2i= (mode 1 and 2 are orthogonal), and the other 
mutual coupling coefficients are identical: hCij = hCji = n 
(i=l,2; j=3). In the coupling region we ignore the dis- 
persion of group velocities and assume that Vg^g^ and 
n are constants. By introducing the following vari- 
ables: ai{z,t) = Ai{z,t)e-'^^, a2{z,t) = A2{z,t)e-'^\ 
as{z^t) = As{z^t)e'^^^ and applying the Fourier trans- 
formation, we obtain the propagation constants of three 
eigenmodes: ki^2 = (a ± iy^—Svg^Vgtz'^ — /3^) /2vgVg^, 
ks = {oj -\- S)/vg where a = Vg{uj — S) -\- Vgs{oJ + S) and 
P = Vg{uj - 5) - Vgs{(jj + 5). When —Svg^VgHi'^ > /3^, ki^2 
is a conjugate pair, indicating the existence of a spec- 



tral gap, while ks corresponds to eignemode a{k,u) = 
di{k^uj) +a2(/c,cc;), where di{k^uj) and a2{k^uo) denote 
orthogonal circularly polarized modes. Thus, k^ corre- 
sponds to linearly polarized HEn mode, which is not 
coupled to the TMqi mode. This mode makes the gap 
polarization dependent. Fig. [TJb) shows the results ob- 
tained using TCMT of S = when values of Vg^gs are 
taken from Fig. [2l^a). 
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FIG. 3. (Color online) (a) Dispersion and (b) losses (imagi- 
nary part of rieff — kz/ko — kzc/uj) of the three eignemodes 
of the two- waveguide structure. Dashed black (7=0), solid red 
and green (7=0.002) curves correspond to modes of conjugate 
propagation constants. Blue curve (dashed curve almost over- 
laps the solid curve as the loss of this mode is comparatively 
low as shown in (b)) correspond to the HEn mode that is 
not coupled to the TMqi mode as shown in (c). (c)-(e) Sz of 
modes at the points (A)-(C) marked in (a), respectively, (f) 
Gap region vs distance between waveguides when 7=0. 

Full numerical simulation results using COMSOL (see 
Fig. [3]) qualitatively agree with TCMT. In the lossless 
case 7 = 0, the spectral gap is defined by a pair of 
complex conjugated propagation constants [see Figs.[3l^a) 
and (b)]. The gap width increases with decreasing the 
distance D [see Fig. [3l^f)], because the coupling coeffi- 
cient becomes larger. When we incorporate some losses 
(7 = 0.002), all modes become complex, and the def- 
inition of width of the gap depends on how far it is 
from the observing point to the source. However, the 
gap width of lossless metal (7 = 0) may still serve as a 
guide and effective approximation as shown in Figs. [3l^a) 
and (b). In addition to the modes of conjugate prop- 
agation constants, there exists one more HEn-like de- 
coupled mode. The energy flow of this mode is mostly 
confined inside £1=4 waveguide [see Fig. [3l^c)]. Thus, 
the gap of the two coupled waveguides is incomplete and 
polarization dependent. To make modes of different pre- 
ferred polarization directions degenerate and obtain a 
full gap, symmetric structures could be used^— One 
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FIG. 4. (Color online) (a) Schematic of the four- waveguide 
structure with Czv symmetry. The distance between si,2 
waveguide to £3,4,5 waveguides is D. (b) Dispersion and (c) 
losses of three eignemodes. Blue curve corresponds to the 
E3 mode. The complete gap region in lossless case is colored 
yellow, (d) Sz of E3 mode at point (A) marked in (b). 



of the options is to utilize four- waveguide C^v structure 
[see Fig. lU^a)]. We use subscripts n = 1,2 to denote 
two HEii modes of circular polarizations and n = 3, 4, 5 
for three TMqi modes. Based on the symmetry and en- 
ergy conservation law in the lossless case, the following 
relations are satisfied for mutual coupling coefficients: 

1^12 = 1^12 = 0, K.im = l^^i = l^lm = ^m2 = K.^ei''^'^-^^' 



for 771, n = 3, 4, 5 and m ^ n. Due to the C^y symmetry, 
eigenmodes a(/c, uo) = X]m=i Pm^imik^ uo) of preferred x 
polarization (/3i = j32) and those of preferred y polar- 
ization {j3i =—^2) should be degenerate. ^^^^^ Thus using 
TCMT we can find five eigenmodes of three frequencies: 
^1,2 = (a ± i^—12vg^VgK?' — /3^) /2vgVg^ (corrcspoud" 
ing to two degenerate pairs of modes), uo = uo^ik) + 2/^2, 
where a = Vgiuj — 5 ^ ^2) + Vg2,{uj -^ 5)^ j3 = Vgiuo — 5 ^ 
^2) — '^g?>{^ + ^) and a;3(/c) is the dispersion of individual 
TMqi mode. Again ki^2 can be a conjugate pair, indicat- 
ing the existence of a gap. uo = uj^{k)^2i<i2 corresponds to 
eignemode (E3 mode) aik^uo) = Xl^=3 ^m(/^,ct;), which 
is a symmetric combination of TMqi modes. The cut off 
frequency of E3 mode is shifted by 2/^2 compared with 
individual TMqi mode. Numerical results from COM- 
SOL are shown in Fig. HJ This allows us to conclude 
that the spectral gap indicated by the yellow region of 
four coupled waveguides becomes polarization indepen- 
dent when E3 mode is cutoff. For larger losses (metal in 
deep ultraviolet regime) the spectral gap still exists, but 
the effective width becomes smaller and eventually dis- 
appears as increasing losses make the differences between 
gap and non-gap region smaller. To enable coupling at 
longer- wavelength regime, where losses of metal is lower, 
one could use dielectric waveguides with higher permit- 
tivities (GaAs for example). 



In summary, we have studied a coupler based on two 
dielectric waveguides in metal involving the coupling of 
backward and forward waves. By using the temporal 
coupled-mode theory we have predicted a spectral gap in 
such a system without a periodic structure. This result 
has been verified by direct numerical simulations. More- 
over, we have demonstrated that a complete polarization 
independent gap can be achieved by using four coupled 
waveguides with C^y symmetry. Similar coupling be- 
tween surface plasmon polaritons (SPPs) can happen in 
metallic-wire structures when the radius is small enough 
to support backward SPPs.^^ However, high losses of 
backward SPPs on metallic wires prevent them from re- 
alistic realizations. We anticipate that by incorporating 
materials with gain and/or nonlineararities, the proposed 
structure can be considered as a new platform for the 
study of gap solitons, optical bist ability, high-Q cavities, 
plasmonic nanolaser in various systems without period- 
icity. 
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